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1. Introduction. 

The aim of this study is to provide a numerical and asymptotic 

description of the structure of planar laminar flames, propagating in a 
medium containing a uniform cloud of fuel-particles premixed with air. 
Attention is restricted here to systems where the fuel-particles first vaporize 
to form a known gaseous fuel, which is then oxidized in the gas-phase. This 
program is supported for the period September 14, 1991 to September 13, 1992. 
Some results of the study is shown in Ref. 1. The work summarized in Ref. 1 
was initiated prior to September 14, 1991 and was completed on February 1992. 
Research performed in addition to that described in Ref. 1 in collaboration 

with Professor A. Linan, is summarized here. 

The model developed here is built on previous asymptotic analysis 
structure of premixed particle-cloud flames [1]. In this previous analysis the 
vaporization rate of the fuel-particles was presumed to be proportional to 
T s l-33 t where T s is the temperature of the fuel-particles. Hence, the activation 
temperature in an Arrhenius-type rate law for the vaporization process is not 
large. For simplicity, the temperature of the fuel-particles was presumed to be 

equal to the gas temperature. The analysis was performed for values of <{>u > 1.0, 
where <t> u is the equivalence ratio based on the gaseous fuel available in the 
particles. For given values of <j> u the analysis yields results for the burning 

velocity and <{>g, where <j>g is the effective equivalence ratio in the reaction 

zone [1]. The analysis shows that even though <j»u > 1.0, for certain cases the 
calculated value of <|>g is less than unity. 

A two-temperature model is described here, where the temperature of 
the fuel-particles is different from the temperature of the gas. However the 
size of the fuel-particle is assumed to be small enough so that its velocity is 
nearly equal to the gas velocity. A general treatment of flame propagation 
supported by volatile fuel-particles should consider both radiative and 
molecular transport mechanisms, the temperature difference and the 
differences in velocity between the gas and the particulates, the 
nonadiabaticity of combustible systems of finite size, as well as the detailed 
kinetics of oxidation and the kinetics of vaporization pyrolysis. This work has 
the somewhat limited aim of examining the interplay of vaporization kinetic 
and oxidative kinetic processes for cases where nonadiabatic and radiative 
transport mechanisms are not substantial. 
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2. Formulation. 

A model is developed to describe steady, one-dimensional, P 1 "*" fl *“ c 

propagation in a combustible mixture consisting of uniformly distributed fuel 

[Zcfes i„ air. The initial number density of the partteles, n u (number of 

particles per unit volume) and the initial radius a u are presumed to be known. 
All external forces including gravitational forces effects are assumed to be 
negligible. The kinetics of vaporization of a fuel-particle is represented by 

the expression 


m = 4ita2Aexp(-^) = -|7t^lp 


0) 


where rh represents the mass of gaseous fuel vaporized per unit time from the 
fuel-particle, a and p s are the instantaneous radius and density respectively o 
a fuel- particle. The rate parameters A and T av which respectively represent 
the frequency factor and the activation temperature of the vaporization 
process Te pLum ed to be known. The ratio T av /T s is presumed to be a large 

quantity." ~ ..... ... ...... V. .... 

Consider a vaporizing fuel-particle in an ambient atmosphere The 
balance equation for the mass fraction of the gaseous fuel vaporizing from the 
fuel-particle Y F , and the energy balance equation can be written as, 
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where r is the radial coordinate. For the gas mixture k is the thermal 
conductivity, c p the heat capacity, and T the temperature. The quantity L F 
k/(pcnD f) is the Lewis number of the gaseous fuel which is presumed to e a 
constant and D F is the diffusion coefficient of the gaseous fuel. Boundary 
conditions for Eqs. 2 and 3 far from the surface of the fuel -particle, and at the 
surface of the fuel-particle can be written as, 

T = T„, Y F = Y P oo as r->o«, (4) 

T = T,.ri,(l-Y F ) = -^^a,r = , 

The values of T ~ Yp >~ will be determined from the solution of the 
equations governing the structure of the laminar flame The temperature o 

the fuel-particle T s can be calculated from the energy balance equation 
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t= a 


( 5 ) 
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where t represents the time, c s represents the heat capacity of fuel-particle, 
and L v represents the heat required to vaporize unit mass of vapor from a fuel- 
particle. The values of c p and c s are presumed to be constant. If the value of k 
is also presumed to be constant, then from integration of Eqs. 2 and 3, using Eq. 
4 to determine the constants of integration, the profiles of Yp and T can be 
expressed as 


and 


Y F = 1 -(1 - Y F „)exp 
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Introducing the definition, 


X = 


U K k r/ r = a’ 
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from Eqs. 5 and 7 it follows that 
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The term on the left side of Eq. 9 denotes the rate of heat transfer from 
the gas to the fuel-particle, with the effective Nusselt's number denoted by 

X/[exp (A.) - 1]. The first term following the equality sign on the right side of 
Eq. 9 denotes the latent heat required to vaporize the fuel-particle, and the 
second term denotes the sensible heat required to raise the temperature of the 
fuel-particle. It follows from Eqs. 1, and 8 that 

X, = ^ A exp . (10) 

Let T v denote a characteristic value of the temperature of the fuel- 
particle for which the value of X evaluated from Eq. 10 is of order unity. For T s 
< T v , the vaporization rate and consequently the value of X, is exponentially 
small. In the limit of small values of X, X/[exp(X) - 1] = 1. Hence the first term 
on the right side of Eq. 9 can be neglected in calculating the value of T s . For T s 
= T v , X is of order unity. If the value of T s is allowed to increase beyond T v , 
then the value of X would become exponentially large and the term on the left 
side of Eq. 9 would become exponentially small, consequently the rate of heat 
transfer from the gas-phase to the solid, will be small. Therefore, for values of 
T s beyond T v , the value of the derivative dT s /dt will be negligibly small, so that 

the value of X remains of order unity. For simplicity, it is presumed that rh = 0 

V 
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for T s < T v and the value of T s does not increase after it attains the value T v . 

The governing equations for the structure of the flame can be written 
as. 


V = p u U u + 


I It ps n u u u (a2 - a 3 ). 
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d(VYoJ k„ d 2 Yoj _ ... Pu 

dx Cp Loj dx 2 ^ P ’ 

(13) 

n “P“ U “ qp+Q(t w F ^ 
dx dx 2 v P 

(14) 


Here the subscript u denotes conditions in the ambient reactant stream. 


The subscript 00 has been dropped from the quantities Yp, Y 02 a °d T appearing 
in Eqs. 12-14. The quantities wp and wo 2 denote the rate of consumption of 
gaseous fuel and oxygen per unit volume, and Qr is the heat released per unit 
mass of gaseous fuel consumed. The independent variable x is related to the 
spatial coordinate x' as. 



(15) 


In the regions where T s < T v , the vaporization rates are assumed to be 
negligibly small hence the radius of the fuel-particle does not change and its 
temperature T s and the quantity q p can be calculated from the relations. 


and 
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When the value of the temperature of the fuel-particle reaches Ty, they 
begin to vaporize and the increase in their temperature is negligibly small. 
The radius of the fuel-particle and the quantity q p can be calculated .(rpro the 
expressions, 


a da _ P“ 
dx Ty Cp Ps V 


T-T v 
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and. 
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qp = m(c p -L v ),m = -p s |x^^. 


(18) 

(19) 
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In deducing the system of equations 12-19 the product pk was presumed 
to be a constant. The chemical reaction between the gaseous fuel and the 
oxidizer is presumed to occur by a one-step process. The system of equations 
12-19 were solved numerically. 

3. Results 


For purpose of illustration calculations were performed assuming that 
the gaseous fuel that evolves from the fuel-particle is methane. The chemical 
kinetic rate parameters characterizing the gas-phase oxidation were chosen 
such that the burning velocity for stoichiometric mixture of gaseous methane 
and air is 37 cm/s. The values of some of the parameters used in the 

calculations arc k u = 0.00035 cal/(cm s K), p u = 0.001135 g/cm^, p s = 1.0 g/cm^. 
T u = 300 K, T v = 600 K, L v /Qr = 0.01. Results of some of these calculations for 
values of a u equal to 0.0007 cm, 0.0010 cm, and 0.0015 cm are shown in Figs. 1 
and 2. Figure 1 shows the the burning velocity u u , as a function of the 

equivalence ratio <}> u , calculated based on the gaseous fuel available in the fuel- 
particles. Figure 1 shows that steady flame propagation is possible in fuel-rich 
combustible mixtures. In Fig. 2 the structure of the flame is plotted at a value 

of <|> u = 1.89 and a u = 0.0010 cm. The various non-dimensional quantities 
appearing in this figure are defined as. 


p u ti u c p T-T u fl _ T S -T„ 

I ^ ^ rp rp 9 S pp rp 9 

Ab “ 1 u lb“*u 

a ^ yF= ^’ yo>= W 


( 20 ) 


where Tb is the adiabatic flame temperature. Interpretations of results similar 
to those shown in Fig. 1 and 2 are currently in progress. 

4. Future Research. 


Future research is directed toward performing numerical calculations 
using the system of equations 12-18 over a wide parametric range. Asymptotic 
description of the flame structure will also be attempted. 
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Figure 1. The burning velocity u y (cm/s) as a function of the equivalence 
ratio <t> u , for various values for the initial radius of the fuel-partical a u - 



Figure 2. The structure of the flame for <|) u = 1.89 and a u = 0.0010 cm. 
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